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ABSTRACT 

Intermediate resolution (0.5-1 A) optical spectroscopy of the cataclysmic variable SY 
Cnc reveals the spectrum of the donor star. Our data enable us to resolve the orbital 
motion of the donor and provide a new orbital solution, binary mass ratio and spectral 
classification. We find that the donor star has spectral type G8±2 V and orbits the 
white dwarf with P = 0.3823753 ± 0.0000003 day, K 2 = 88.0 ± 2.9 km s^ 1 and 
V smi = 75.5 ± 6.5 km s . Our values are significantly different from previous works 
and lead to q = Mi/M\ = 1.18 ± 0.14. This is one of the highest mass ratios known 
in a CV and is very robust because it is based on resolving the rotational broadening 
over a large number of metallic absorption lines. The donor could be a slightly evolved 
main-sequence or descendant from a massive star which underwent an episode of 
thermal-timescale mass transfer. 
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1 INTRODUCTION 

Mass transfer in cataclysmic variables (CVs hereafter) is 
driven by angular momentum loss mechanisms. Magnetic 
braking is considered the main driver for orbital periods 
longer than ~3 h, although its actua l prescription is un- 
der debate (e.g. Ilvanova fe Taamll2003h . As a result of mass 
transfer, the radius of the donor star is adjusted (1) on a 
thermal time-scale to try preserving thermal equilibrium and 
(2) on a dynamical time-scale to return to hydrostatic equi- 
librium. The competition between the changing stellar and 
Roche lobe radii sets the mass transfer rate and hence its 
stability. For large mass ratios the binary could be unstable 
agains t rapid ma ss transfer and a common envelope may 
ensue (|Kindll988l ). 

SY Cnc is one of the brightest Z Cam-type d warf no- 
vae, e xhibiting regular outbursts every ~27 days. IShafterl 
(|1985P reported an orbital period P = 0.380 ± 0.001 d and 
q — M2/M-1 — 1.13±0.35 which was interpreted as evidence 
for unstable mass transfer. Unfortunately, q was derived in- 
directly using observed properties of the Ha emission line 
and this method is not entirely reliable because the accre- 
tion disc's emissivity is often found to be non-axisymmetric 
and contaminated by hot-spots. The best way to constrain 
q is by combining the radial velocity semiamplitude of the 
donor star K2 w ith the rotational bro adening of its absorp- 
tion lines V sin i ijWade fe Hornlll988l V The issue of the high 
q value was revisited by Smith et alj ([2005) who measured 
K2 = 127 ± 23 km s _1 from a skew map of the Na I IR 



doublet. However, the poor spectral resolution hindered a 
direct determination of V sin i and hence K2 was combi ned 
with K\ (from the wings of the Ha line, IShaftenll983h to 
yield a revised q — 0.68 ± 0.14. They also classify the donor 
star as M0, which must be substantially evolved in order to 
fill its 9.1 h Roche lobe. 

Here we present new higher resolution spectra of SY 
Cnc demonstrating that the spectrum of the donor star is 
clearly detected at optical wavelengths. We update the or- 
bital parameters and present the first determination of the 
rotational broadening of the donor star, which removes the 
lingering uncertainty on the true mass ratio. We also find 
that the spectrum of the donor is best fitted by a late G 
star and discuss the observed properties in the light of mass- 
transfer stability and the evolution of SY Cnc. 



2 OBSERVATION AND DATA REDUCTION 

SY Cnc was observed over several nights between 1992 and 
2008 using the Intermediate Dispersion Spectrograph (IDS) 
attached to the 2.5 m Isaac Newton Telescope (INT) at the 
Observatorio del Roque de Los Muchachos. The 27 spectra 
were always obtained in the Ha region but using three grat- 
ings and slightly different instrumental settings (slit widths 
and central wavelengths), which resulted in spectral resolu- 
tions in the range 25-70 km s" 1 . The nights were allocated 
to various scientific programs and hence only a few spectra 
of SY Cnc were obtained per night. A log of the observations 
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Table 1. Log of the observations. Il976f ) and hence we can conclude that 0.3824 day is by far 

the most significant period. A least-squares sine-wave fit to 
the radial velocities, using P — 0.3824 day as input param- 
eter, yields the following parameters 



Date 


Wav. Range 


Exp. time 


Dispersion 


Resolution 




A A 


(s) 


(A pi*- 1 ) 


(km s _1 ) 



10/03/1992 


6025-6975 


2400 


0.79 


70 


13/02/1998 


6290-6800 


600,900 


0.53 


411 


1-3/05/1999 


6290-6800 


2x600,6x900 


0.53 


49 


22-27/04/2003 


6100-6780 


10x600 


0.30 


25 


20-22/10/2008 


5575-7139 


3x300,3x600 


0.64 


53 



7 = -0.2 ± 2.5 km s" 1 
P = 0.3823753 ± 0.0000003 d 



is presented in Table [T] Spectra of the stellar templates 61 
Cyg A & B (K5V and K7V) were also observed during the 
2008 run for the purpose of radial velocities and rotational 
broadening analysis. In addition, eight templates of spectral 
types G6-M0V were collected during the 2003 campaign and 
another K5V in 1999. 

The images were bias corrected and flat-fielded, and the 
spectra subsequently extracted using conventional optimal 
extraction techniqu es in order t o optimize the signal-to-noise 
ratio of the output (|Hornell986l ). Every target was bracketed 
with observations of a comparison CuAr+CuNe arc lamp 
and the pixel-to-wavelength scale was derived through poly- 
nomial fits to a large number of identified reference lines. 
The final rms scatter of the fit was always <l/30 of the 
spectral dispersion. 



3 THE ORBITAL SOLUTION 

We rectified the 27 individual spectra by subtracting a low- 
order spline fit to the continuum, after masking out the main 
emission and atmospheric absorption lines. The spectra were 
subsequently rebinned into a uniform velocity scale of 36 km 
s _1 pix -1 . The K5V template 61 Cyg A was broadened to 78 
km s _1 to match the width of the donor photospheric lines 
(see Sect. 4). Every spectrum of SY Cnc was then cross- 
correlated against the broadened template in the spectral 
regions free from emission and telluric absorption features. 
Radial velocities wer e extracted following the method of 
iTonrv fc Davis! (|l979f ). where parabolic fits were performed 
to the peak of the cross-correlation functions, and the un- 
certainties are purely statistical. Since the orbital period is 
poorly constrained to 0.380 ±0.001 d we performed a power 
spectrum analysis on the radial velocities in the range 0.1-2 
days, and the results are displayed in Fig. [T] Here we have 
rescaled the errorbars by a factor 1.3 so that the minimum 
xt is 1.0. The periodogram is dominated by strong aliasing 
due to the sparse sampling of our observations. In order to 
test the significance of the different peaks above the noise 
level we performed a Monte Carlo simulation. Synthetic % 2 
periodograms were computed from a large population (10 ) 
of velocities randomly picked from a white-noise distribution 
and with identical time sampling as our data. A 4-cr signifi- 
cance level is defined by the 99.99 per cent of the computed 
xt values and this is indicated in Fig. [J by a horizontal 
dashed line. Most of the peaks between frequencies 2.5-2.7 
day -1 are under the line and hence periods in the range 
0.37-0.40 days are significantly above noise at the 99.99 per 
cent level. The two deeper peaks correspond to 0.3824 day 
and 0.3837 day and have xt — 1-0 an d 2.3 respectively for 
24 degrees of freedom. A 4-<r significance l evel around the 
minimm peak will exclude the second peak l|Lampton et al.l 



To = 2451300.254 ± 0.002 

K 2 = 88.0 ± 2.9 km s" 1 

where To corresponds to the Heliocentric Julian date of the 
inferior conjunction of the donor star. All quoted errors are 
68 per cent confidence. The systemic velocity 7 has been 
corrected from the rad ial velocity of 61 Cyg A, that we take 
as -64.3 ±0.9 km s" 1 l| Wilson II 1953 ). Note that a sinewave 
fit fixing P = 0.3837 day also yields K 2 = 87.3 ± 3.3 km 
s _1 , an indication that our K2 value is robust and its error 
realistic, irrespectively of the true value of the orbital period. 
The same is found when other peaks around the minimum 
are taken. 

Fig. [5] displays the radial velocity points folded on our 
favoured orbital period P = 0.3823753 day together with 
the best sine fit solution. Our K2 ve locity disagrees with the 
one reported bv lSmith et al. I l|2005h which was obtained us- 
ing the 8190-A Na 1 doublet. We note that the Na 1 doublet 
can be quenched b y heating effects, as opposed to the metal- 
lic lines used by us l|Martin et al.|[l989l) . If this were the case, 
the light centre of the Na I lines would be displaced towards 
the back side of the star. The effects of irradi ation can be es- 
timate d using the K-correction approach of Wa de fc Hornel 
l| 19881 ) AK2/K2 = /(r 2 /a)(l ± q), where AK 2 is the incre- 
ment in K2 velocity due to irradiation, T2 the radius of the 
donor star, a the binary separation and / the fractional dis- 
placement of the absorption line's site with respect to the 
center of mass of the star. In the extreme case, when the 
Na 1 absorption is completely supressed from the irradiated 
hemisphe re, / = 4/(37r). Thus, replacing T2/a by Eggleton's 
equation (Eggleton 1983) and adopting q — 1.18 (see next 
section) we find AK2 = 32 km s . This difference is just 
enough to accommodate Smith et al's value with ours, for 
maximum quenching of the Na I lines. Therefore, we con- 
clude that Smith et al's K2 is likely a significant overesti- 
mate due to irradiation while our value is much less affected 
by these effects. 



4 THE ROTATIONAL BROADENING AND 
SPECTRAL CLASSIFICATION 

In order to measure the rotational broadening of the donor's 
absorption features we have focused on the 10 highest resolu- 
tion spectra gathered during the 2003 campaign. We broad- 
ened the G6-M0 templates from 5 to 100 km s _1 in steps 
of 5 km s , using a Gray profile (|Gravlll992l ) and contin- 
uum limb-darkening coefficients appropriate for every spec- 
tral type and our wavelength range. The broadened tem- 
plates were multiplied by factors / < 1, to account for the 
fractional contribution to the total light. These were sub- 
sequently subtracted from the Doppler corrected average of 
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donor star around the minimum value. The horizontal dashed 
line indicates the 99.99 per cent significance level above noise. 
The deepest peak corresponds to P = 0.3823753 d. 
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Figure 2. Radial velocity curve folded on P = 0.3823753 day. 
The best sine-wave solution is overplottcd. 



SY Cnc, obtained using our orbital solution above. The op- 
timal broadening, based on a \ 2 test on the residuals, is 
always found in the range 77-80 km s _1 . The mean of these 
determinations is 78.4 ± 3.5 km s _1 . The quoted error cor- 
responds to 1-ct and has been estimated by fitting a cubic 
function to the x 2 vs V sin i curve and searching for Xmm + !• 
A potential source of systematics in this calculation is 
the assumption of continuum limb-darkening coefficient be- 
cause absorption lines in late-type stars are expected to have 
smaller core limb-darke ning coefficients than the continuum 
l|Collins fc Truaxlll995l ). We have tested this by repeating 
the above analysis using zero limb-darkening as a secure 
lower limit and obtain Vsini = 72 ± 3 km s -1 . Therefore, 
we decided to adopt 75.5 ± 6.5 km s _1 as a conservative 
value for V sin i hereafter. The minimum x 2 a lso constrains 
the donor's spectral type to G9-K3V. The same analysis 
was independently performed using the lower resolution data 



1 i , 1 1 




1 






# K7 H 




K3 

• 




GO 


K4 

• 




G1.5 

K1.5 

G6.5 K0 
• • 




i 


3.75 3.7 


3.65 


3.6 


log T eff 


(K) 





Figure 3. Spectral type classification of the donor star. The min- 
imum x 2 is well constrained between spectral types G6.5V and 
K0V. 



from 1999 and 2008, which yield consistent results at 80 ± 5 
and 77 ± 5 km s~ , respectively. In the latter case the K5V 
template produces a lower minimum \ 2 than the K7V, an 
indication that the donor star is earlier than K5. 

To further constrain the spectral type we employed 
a set of high resolution (R=A/AA=55000-11000) spectra 
covering G0-K7V from Ecuvillon et al. (2004, 2006) and 
iMontes fc Martini (|l998T l. These are HD 39091 (GOV), HD 
30495 (G1.5V), HD 43162 (G6.5V), HD69830 (K0V), HD 
17925 (K1.5V), HD 222237 (K3V), HD 216803 (K4V) and 
HD 157881 (K7V). These stars have very accurate stellar 
parameters logg,T c ff and metallicities so they can provide a 
more robust spectral type classification. These were down- 
graded to the instrumental resolution of the 2008 INT data 
and broadened using our Vsini estimate. Fig. [3] displays 
the x 2 °f the optimal subtraction as a function of the tem- 
plate's spectral type. The minimum is found for G6.5V and 
K0V and thus we classify the donor star as a G8 ±2V. 
This is again remar kably different from the M0 suggested by 
ISmith et all Jg005) wherea s it is in better agreement with 
the work of lHarrison et afl l|2004l ) who propose an early G 
star based on IR spectra. Fig. [4] presents the 2008 average 
spectrum of SY Cnc, Doppler corrected in the rest frame 
of the donor star, together with one of the best templates. 
The latter must be scaled by a factor 0.42T0.02 to match 
the depth of the absorption features in SY Cnc and hence 
we conclude that the donor star contributes 42 per cent to 
the total flux in the Ha region. 
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Figure 4. Doppler corrected average of SY Cnc in the rest frame 
of the donor star together with the G6.5V template and the resid- 
ual after an optimal subtraction. The template has been broad- 
ened to match the Vsini of the donor star. 



5 DISCUSSION 

Since the donor star is filling its Roche lobe and synchro- 
nized, we can combine our determination of K2 and V sin i 
to constrain the binary mass ratio through the equation: 



Vsvsxi = K 2 (1 + q) 



0.49 q 



2/3 



0.6 q 2 / 3 +ln(l + g 1 /3 



(1) 



This uses Egglet on's expression for the effective radius of 
the Roche lobe l|Eggletonl [m$L_ which is more accurate 
than Pacyriski's approximation (|paczvriskilll97ll ). Substitut- 
ing our values of K2 and Vsini into equation (1) we find 
q — 1.18 ± 0.14, where the error has been estimated using a 
Monte Carlo simulation. The input parameters Vsini and 
K 2 have been picked up randomly from normal distributions 
whose mean and sigma correspond to the observed values 
and the process is repeated 10 5 times. This is one of the high- 
est, accurately measured, mass ratio i n a CV, only compara- 
ble to V363 Aur with q = 1.17±0.07 l|Thoroughgood et al. I 
120041 ). For instance, the 7.10 e dition of the Ritter & Kolb 
catalogue (|Ritter fc Kolbl 2003) lists seven systems in 731 
entries with q > 1. However, only V363 Aur and AC Cnc 
possess robust q values because these are derived through a 
measurement of K2 and V smi. The remaining five should 
be regarded as tentative given the difficulty in deriving stel- 
lar masses in the absence of donor stellar features. 

Theoretical studies of binary evolution show that stable 
mass transfer depends on both the thermal and dynamical 
response of the donor star to mass loss. This is described by 
the thermal and adiabatic radius-mas s exponents, co mputed 
for the case of ZAMS donor stars re.g. |Politanolll996h . When 
compared with the tidal radius-mass exponent (i.e. the log- 
arithmic derivative of the Roche lobe radius with respect to 
the donor mass) they define a critical mass ratio for stable 
mass transfer. The radius-mass exponents, in turn, depend 
on the donor's mass so a combination of (q,M2) values is 
required to fulfil the condition for stable mass transfer. 

SY Cnc presents a remarkable case because of its high 
(and accurate) mass ratio which places it close to the ab- 
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Figure 5. Binary mass ratio versus donor mass. The solid and 
dashed lines represent the critical mass ratio for adiabatic (dy- 
namical) and thermal mass-transfer stability, respectively, for 
ZAMS donor stars. The thick horizontal lines mark our limits 
to q. 



solute upper li mit of mass transfer stability computed by 
IPolitanoT (|l996t ). Fig. [5] shows the critical mass ratio curves 
for dynamical (solid line) and thermal stability (dashed line) 
as a fu nction of the donor m a ss, up dated after figs. 1 and 
2 from ISmith fc Vande Puttel (|2006l ). For SY Cnc to lie in 
the stable mass-transfer region (i.e. below the two curves) 
our central value q = 1.18 yields M 2 ~0. 77-1. 54 M . This 
would imply a white dwarf mass Mi ~ 0.65 — 1.31 Mq and, 
since Mi sin 3 i = (1 + q) 2 PK'i /2ttG = 0.13 ± 0.02 Mq, then 
i ~ 26 — 38°. However, the theoretical stability limits are 
computed under the assumption of ZAMS donor stars which 
is not necessarily true. Donor stars in CVs have on aver- 
age later spectral types than field m ain-sequence stars , an 
indicati on of radius expa n sion ( e.g. iKniggel l2006t ) . In par- 
ticular, iBeuermann et all (Il998l ) showed that donor stars 
with P or b > 5 — 6 h were likely nuclear-evolved at the start 
of mass transfer. This alters the helium-to-hydrogen com- 
position and hence the response of the star to mass loss. 
Therefore, the critical mass ratios shown in Fig. [5] should be 
regarded as mere indicative limits. 

Our determination of the spectral type also provides a 
new ingredient to constrain the evolution of the donor star 
in SY Cnc. A G8 ZAMS has M 2 = 0.81 M Q which would 
place it reasonably well within the stability region. A donor 
with such mass would fill ~ 87 per cent of its 9.1 h Roche 
lobe. The conclusion seems to be that the companion star is 
a slightly evolved main sequ ence, in contrast with previous 
claims bv lSmith et all l|2005h . 

Alternatively, SY Cnc could be descending from a pro- 
genitor binary with a more extreme mass ratio which under- 
went a phase of thermal-timescale mass transfer. The impor- 
tance of this ev olutionary path has been addressed by sev- 
eral papers (e.g.lSchenkerfc King 2002; Podsiadlo wski et al.1 
120031 ; iKolb fc Willemsl 120051 ) and successfully explains the 
cooler-than-main-sequence donor stars typically observed at 
> 5 — 6 h orbital periods. Natural consequences of this sce- 
nario are larger white dwarf masses and evi dence for CNO 
processed material, both seen in AE Aqr (Sch enker et al.l 
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2002). The detection of nuclear processed material in UV 
spectra of SY Cnc would clarify the evolutionary history of 
SY Cnc and whether it is the outcome of a thermal mass- 
transfer episode. Further, more higher resolution observa- 
tions will provide better constraints on V sini, K2, the stel- 
lar masses and the evolution of the donor star in this impor- 
tant CV. 



6 ACKNOWLEDGMENTS 

We thank the anonymous referee for the very helpful 
comments which improved the quality of the manuscript. 
MOLLY software developed by T. R. Marsh is gratefully 
acknowledged. Partly funded by the Spanish MEC under 
the Consolider-Ingenio 2010 Program grant CSD2006-00070: 
first science with the GTC. Based on observations made with 
the INT operated on the island of La Palma by the Isaac 
Newton Group in the Spanish Observatorio del Roque de 
Los Muchachos of the Instituto de Astroffsica de Canarias 
(I AC). 



Smith R.C., Vande Putte D., 2006, The Observatory, 126, 
38 

Thoroughgood T.D., Dhillon V.S., Watson OA., Buckley 
D.A.H., Steeghs D., Stevenson M.J., 2004, MNRAS, 353, 
1135 

Tonry J., Davis M., 1979, AJ, 84,1511 
Wade R.A., Home K., 1988, ApJ, 324, 411 
Wilson R.E., 1953, in General Catalogue of Stellar Radial 
Velocities, Carnegie Inst. Washington DO. Publ. 601, 



REFERENCES 

Beuermann K., Baraffe I., Kolb U., Weichhold M., 1998, 

A&A, 339, 518 
Collins G.W., II, Truax R.J., 1995, ApJ, 439, 860 
Ecuvillon A., Israelian G., Santos N. O, Mayor M., Villar 

V., Bihain G., 2004, A&A, 426, 619 
Ecuvillon A., Israelian G., Santos N. O, Shchukina N. G., 

Mayor M., Rebolo, R., 2006, A&A, 445, 663 
Eggleton P.P., 1983, ApJ, 268, 368 

Gray D.F., 1992, The Observation and Analysis of Stellar 

Photospheres. CUP, Cambridge 
Harrison T.E., Osborne H.L., Howell S.B., 2004, AJ, 127, 

3493 

Home K., 1986, PASP, 98, 609 

Ivanova N., Taam R.E., 2003, ApJ, 599, 516 

King A.R., 1988, QJRAS, 29, 8 

Knigge C, 2006, MNRAS, 373, 484 

Kolb U., Willems B., 2005, The Astrophysics of Cata- 
clysmic Variables and Related Objects, ASP Conf. Series, 
330, 17 

Lampton M., Margon B., Bowyer S., 1976, ApJ, 208, 177 
Martin J.S., Friend M.T., Smith R.C., Jones D.H.P., 1989, 

MNRAS, 240, 519 
Montes D., Martin E.L.., 1998, A&AS, 128, 485 
Paczyhski B., 1971, ARA&A, 9, 183 

Podsiadlowski Ph., Han Z., Rappaport S., 2003, MNRAS, 

340, 1214 
Politano M., 1996, ApJ, 465, 338 
Ritter H., Kolb U., 2003, A&A, 404, 301 
Schenker K., King A.R., 2002, The Physics of Cataclysmic 

Variables and Related Objects, ASP Conf. Series, 261, 242 
Schenker K., King A.R., Kolb U., Wynn G.A., Zhang Z., 

2002, MNRAS, 337, 1105 
Shaffer A.W., 1983, PhD thesis, Univ. California 
Smith DA., Dhillon V.S., 1998, MNRAS, 301, 767 
Smith R.C., Mehes O., Vande Putte D., Hawkins N.A., 

2005, MNRAS, 360, 364 



